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Abstract. The energy-momentum tensor of a supersonic heavy quark jet moving 
through a strongly-coupled J\f = 4 SYM plasma is analyzed in terms of first-order 
Navier-Stokes hydrodynamics. We focus on the near zone "head" region close to the 
heavy quark where deviations from hydrodynamic behavior are expected to be large. 
For realistic quark velocities, v — 0.99, we find that the non-hydrodynamical head 
region is confined at a narrow pancake surrounding the jet with a width ~ 1/T and 
transverse radius ~ 3/T. Outside this jet head region, the jet induced stress is well 
approximated by Navier-Stokes hydrodynamics. 

PACS numbers: 25.75.-q, 11.25.Tq, 13.87.-a 
1. Introduction 

The suppression of highly energetic particles observed at the Relativistic Heavy Ion 
Collider (RHIC) [1] suggests that the matter created in the first moments after the two 
heavy ions have collided is a strongly-coupled deconfined plasma of quarks and gluons [2] . 
The degree of thermalization of the soft degrees freedom seems to be consistent with 
the perfect fluid hypothesis where the strongly-coupled quark-gluon plasma (sQGP) 
has a shear viscosity to entropy density ratio [3] of the order of the string-inspired 
lower bound rj/s = 1/ (47r) [I] derived using the anti-de Sitter /Conformal Field Theory 
(AdS / CFT) correspondence [5] . This correspondence provides a framework to calculate 
nonpertubatively observables associated with strongly-coupled M = 4 supersymmetric 
Yang-Mills (SYM) theories in terms of dual weakly-coupled (super)gravity. 

If the RHIC plasma is indeed strongly-coupled, one has to rely on non-perturbative 
methods, such as AdS/CFT, to describe its properties. At temperatures sufficiently 
above the critical temperature, the sQGP appears to display some of the inherent 
properties of M = 4 SYM plasmas such as conformal invariance. Given the power of 
AdS/CFT to compute nonequilibrium dynamical quantities, it is natural to investigate 
whether jet quenching phenomena can be explained in terms of the dual gravitational 
concepts of AdS/CFT. 

In this paper we therefore study the SYM analogous of jet quenching, i.e., the 
problem of a heavy quark moving through a strongly-coupled Af = 4 SYM plasma. 
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We use the setup proposed in [6] where an infinitely massive heavy quark moving at a 
constant speed v along the comoving X\ — X— vt axis in the SYM plasma corresponds to 
one of the endpoints of a classical string that lives in a 5-dimensional space defined by an 
AdS5-Schwarzschild black (brane) hole geometry. In this steady-state solution, the quark 
has been moving since t — > — oo and can be found at the origin at t = 0. According to the 
AdS/CFT correspondence, the disturbances in the 5-dimensional metric caused by the 
string can be used to obtain the full non-perturbative result for the energy-momentum 
tensor of the system plasma+heavy quark at finite temperature T. In Ref. [7j, the total 
energy-momentum tensor of the system was computed analytically to leading order in 
an expansion in powers of inverse momentum T/K <C 1 and then transformed back to 
coordinate space. This tensor describes the energy and momentum disturbances induced 
in the plasma in a spacetime region close to the heavy-quark in the limit where both 
the t'Hooft coupling constant of the plasma A and the number of colors N c are large. 
Note that the low-energy excitations found in the linearized analysis of Ref. [8], i.e, the 
sound and diffusion modes, only appear in the opposite limit where T/K 1 [9j. In the 
present work, we compare the analytical result for the stress tensor derived in [7] with a 
first-order Navier-Stokes ansatz for the energy-momentum tensor in order to determine 
the distance within which the wake left behind by the heavy quark is thermalized (see 
also [IDJ and pLlJ). 



2. Comparison with first-order Navier Stokes hydrodynamics 

To first order in a gradient expansion, the Navier-Stokes stress tensor is 

= (p + p) U^U U + pg^ + IV (1) 

where p is the local energy density, p is the isotropic pressure, and If 1 = (U°, U) is the 
flow 4- vector (in our metric U fJ- U tl = —1). The dissipative contribution to the tensor 
above is given by the shear tensor (the bulk viscosity of conformal plasmas is identically 
zero) 

W v = -7] ( W + d v U» + U»U a d a U u 

+ U u U a d a U») + \r) A"" (d a U a ) , (2) 

where A» u = gW + U^U V is the local spatial projector and r] = (N c - 1) 2 ttT 3 /8. Note 
that we kept the nonlinear terms in U^ 1 in the definition of IT^. The flow velocity 
associated with the full T^ u of |7J can be obtained by boosting the tensor to its local 
rest frame (Landau frame) (TDJ. The deviation from Navier-Stokes hydrodynamics can 
then be studied by defining the tensor 

Zfiv — — n^j, (3) 

and checking its components in the Landau frame. If (T^j (components evaluated in 
the Landau frame come with the L subscript) can be described by Navier-Stokes, we 
should obtain 

(Z U ) L = (Z 22 ) L = (Z 33 ) L = I (Z 00 ) L , (Z ij ) L = 0. (4) 
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Figure 1. (Color online) is defined in Eq. [3] Ratio values outside of the limits 
were color-coded as the limits. 
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Figure 2. (Color online) Comparison between the full numerical solution for U and 
the leading order result, i.e., \U — S , /(4P )|- 

3. Results 

One can see in Fig. [T] that the magnitude and structure of the discrepancy between 
Navier-Stokes and T^ v remains approximately the same for all the components of (Z >1V ) L . 
We have set N c = 3, A = 3ir, and v = 0.99. Outside the region where —3 < XiirT < 3 
(longitudinal direction) and < X p tcT < 10 (transverse direction), the system can 
be described by the first-order Navier-Stokes energy-momentum tensor to an accuracy 
of roughly 90%. Inside that region, however, the discrepancy quickly increases and 
eventually diverges. In Fig. [2j one can see that in the region where Navier-Stokes 
provides a good description of the system the numerically calculated velocity flow U 
matches the leading order result in l/N c , U = S / (4P ), where Si = —T^ and Po = r] ixT . 
Therefore, in this region the system is more accurately described by a linearized version 
of the first-order Navier-Stokes equations. 
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4. Conclusions 

In this paper we have compared the analytically obtained energy-momentum tensor of 
an infinitely massive quark moving through a strongly coupled M = 4 SYM plasma with 
a non-ideal tensor given by the first-order Navier-Stokes ansatz. We have found that 
the energy deposited by the quark is thermalized in the region defined by |Xi|7rT > 3 
and X p nT > 1, which excludes the nonequilibrium jet head. In fact, the system can 
be described by linearized first-order Navier Stokes hydrodynamics down to distances 
of roughly 3/(ttT). Thus, assuming that a strongly coupled M = 4 SYM plasma can 
capture the main features of the sQGP, we should expect that the energy deposited by 
jets induces almost perfect fluid hydrodynamic response of the sQGP outside the head 
region. 
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